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NUMERICAL METHOD FOR DESIGN OF MINIMUM-DRAG 
SUPERSONIC WING CAMBER WITH CONSTRAINTS ON 
PITCHING MOMENT AND SURFACE DEFORMATION 

By Russell B. Sorrells and David S. Miller 
Langley Research Center 

SUMMARY 

This report presents a numerical method, based on linearized theory, for designing 
minimum -drag supersonic wing camber surfaces of arbitrary planform for a given lift, 
with options for constraining the pitching moment and/or the surface deformation at the 
trailing edge of the root chord and for selecting any desired combination of eight specified 
wing -loading distributions to be employed in the optimization procedure. Two examples 
are presented to illustrate applications of the method. The results indicate that relatively 
small drag penalties are incurred in designing wings to be self -trimming and to have a 
reasonable camber surface. 


INTRODUCTION 

The performance benefits of twist and camber when applied to swept wings with sub- 
sonic leading edges have been demonstrated both theoretically (refs. 1 to 4) and experi- 
mentally (refs. 5 to 7). Frequently, however, the theoretical solution for optimum twist 
and camber results in impractical camber -surface slopes near the wing root. Further- 
more, the resulting wing may or may not be self -trimming at the design lift coefficient. 

The wing-design procedure presented in this report is a logical extension of the 
numerical method based on linear theory described in reference 3 and is intended to pro- 
vide a more versatile system for the design of supersonic wings. Whereas the method 
as applied in reference 3 determines a combination of three specified wing-loading dis- 
tributions that produces minimum drag for a given lift, the present design system provides 
the user the option of selecting any combination of eight specified wing-loading distribu- 
tions for minimizing drag with a given lift, zero pitching moment about a given reference 
point, and a given ordinate at the trailing edge of the root chord. 

To illustrate the application of the present method and to give some indication of 
benefits to be gained, two examples are presented. 



SYMBOLS 


A(L,N) 

Ai 

b 

c 

c 

C D 


C D,ij 

C L 

c L,des 

Cm 

AC p 

k=a 

k=l 

L,N 
L* ,N* 
l 

M 

q 

R 


weighting factor applied to leading-edge grid element (see eq. (3)) 

load strength factor for ith loading 

wing span 

local chord 

mean geometric chord 

drag coefficient 

drag coefficient of interference between ith and jth specified loadings 
(see eq. (20) of ref. 3) 

lift coefficient 

design lift coefficient 

pitching-moment coefficient, Pitchin g moment 

qSc 

lifting pressure coefficient 


k=a i=n j=n . 

X k^k + C D = ^ X k^k + | X Z AiA 3 C D,ij 

k=l i=l j=l 

designation of influencing grid elements (see fig. 3) 
designation of receiving grid elements (see fig. 3) 
length of wing root chord 

Mach number , . 

free-stream dynamic pressure 

average value of influence function within a grid element (see eq. (2)) 
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s 


reference wing area 


x,y,z Cartesian coordinates (see fig. 2) 

x f distance from wing leading edge measured in x-direction 

z c (x,y) camber -surface z-ordinate at x,y 
z r camber -surface z-ordinate at wing-root trailing edge 


0 = 



1 


''-m 

''-z 

<?k 

Subscripts: 


Lagrange multiplier for lift constraint 
Lagrange multiplier for moment constraint 
Lagrange multiplier for z r constraint 
constraining relation 


1,2, 3,. ..,8 

a 

i,j 

k 

le 

max 

min 


specified loadings 

number of constraints 

ith and jth specified loadings 

kth constraint relation 

leading edge 

maximum 

minimum 


n 


number of specified loadings 



NUMERICAL DESIGN METHOD 


Objectives 

The primary objectives of the present design method are: (1) To achieve low trim 
drag through use of a pitching -moment constraint, (2) to provide for more versatility in 
determining an optimum camber surface by increasing the number of specified loadings, 
and (3) to insure practical camber surfaces by imposing a surface -deformation constraint. 
When insufficient pitching moment is generated by the unconstrained design procedure 
(see fig. 1), additional pitching moment must be acquired by deflecting control surfaces. 
When the control surfaces are close coupled, like those shown in figure 1, substantial 
increases in trimmed drag sometimes result. An alternate approach is to include a 
pitching-moment constraint in the wing twist and camber design procedure. This addi- 
tional constraint will result in a drag at the design Cl (see dashed line in fig. 1) which 
in most cases will be less than that resulting from trimming by control-surface deflection. 

General Procedure 

In general, linearized theory is applied to a wing of arbitrary planform and zero 
thickness lying in the xy -plane, and the camber surface is calculated for each of the eight 
specified loadings (to be discussed in detail later). The wing lift and pitching -moment 
coefficients for each of the eight specified loadings, and the wing interference -drag coef- 
ficients between any two loadings, are calculated. Lagrange's method of undetermined 
multipliers is used to minimize drag by determining for each of the eight loadings the 
weighting factor that satisfies the constraints of a given lift, zero pitching moment, and 
a given z c ordinate at the wing-root trailing edge. 

Camber Surface and Force Data for a Specified Loading 

The wing design procedure presented in this report is an extension of the numerical 
method described in reference 3. The numerical method of determining the camber sur- 
face and force data for a specified loading is identical with that of reference 3, and a less 
detailed discussion of the theory and the numerical method used is presented here. 

Figure 2 illustrates a typical wing planform described by a rectangular Cartesian 
coordinate system. Linearized theory for supersonic flow is employed to define the cam- 
ber surface required to support a specified pressure -loading distribution. Numerical 
implementation of the theory utilizes the concept of Mach boxes, in which a grid system 
(see fig. 3) is superimposed on the Cartesian coordinate system. (Fig. 3 is illustrative 
only; in application many more grid elements would be employed.) The symbols L and 
N denote the position of a particular influencing grid element; L* and N* denote the 
position of the element receiving the influence. The region of integration, originally 
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bounded by the wing leading edge and the forward Mach lines from L* ,N*, now consists 
of a set of grid elements approximating that region, as shown by the shaded area in fig- 
ure 3. The wing surface slope at a point represented by L* and N* may now be found 
by summing contributions of all elements within the influencing region by the expression 
(eq. (10) of ref. 3): 

N=N max L=L*-|N*-N| 


9z c 

9x 


(L*,N*) = - | AC p (L*,N*) + iL 



R(L* -L ,N* -N) A(L ,N) AC p (L ,N) 

(1) 


where the brackets in j^x^ designate the whole -number part of the quantity. The func- 
tion R may be thought of as the average value of an influence function within a grid ele- 
ment at L,N relating the local loading at element L,N to its influence in determining 
the necessary slope at downstream element L ,N . This function R is given by the 
expression (eq. (8) of ref. 3): 


R(L*-L N*-N) = \/(L* - L + 0.5) 2 - (N* - N - 0.5) 1 
(L* - L + 0.5)(N* - N - 0.5) 


_ \/(L* - L + 0.5) 2 - (N* - N + 0.5) 2 (2) 

(L* - L + 0.5)(N* - N + 0.5) 

A graphical representation of this factor is shown in figure 4. Note the small variations 
of the factor in the x- or L-direction in contrast to the drastic variations in the Y- or 
N -direction. 

The term A(L,N) in equation (1) is a weighting factor which eliminates the neces- 
sity of accepting or rejecting complete block elements and thus permits a better defini- 
tion of the wing leading-edge shape, which has been found to be extremely critical. The 
factor A(L,N) takes on values from 0 to 1 given by 

A(L,N) = 0 ( L - x le = °f 

A(L,N) = L-x le (0 < L - x le < 1) / (3) 

A(L,N) = 1 ( L - x le = *) 

The section lift and pitching -moment coefficients for each of the eight specified 
loadings and the section interference -drag coefficients between any two loadings i,j are 
then calculated. The wing lift coefficient Cl^ and pitching -moment coefficient C m j 
for each of the eight specified loadings and the wing interference -drag coefficients Cj) jj 
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between any two loadings i,j are calculated by spanwise summations of the section data. 
(See eqs. (16), (17), (18), and (20) of ref. 3.) 

The Eight Specified Wing-Loading Distributions 

The wing design method described in reference 3 employs the following three speci- 
fied wing-loading distributions: uniform, linear chordwise, and linear spanwise. In the 
present method, five additional wing-loading distributions are provided so the drag mini- 
mization procedure will have more versatility in computing the optimum combination of 
loadings. This is particularly important in the present method because of the more 
stringent requirement of satisfying three constraints rather than just one. The eight 
specified wing -loading distributions presently available are illustrated in figure 5 for a 
delta wing, and the advantages and disadvantages of applying these loading distributions 
are discussed in the two illustrative examples. 

Optimum Combination of Loadings 

In reference 3 Lagrange’s method of undetermined multipliers has been applied to 
the problem of selecting the combination of three loadings that yields a minimum drag 
for a given lift. In the present report this technique is again employed to minimize drag 
by selecting the combination of eight loadings which satisfies constraints placed on lift, 
pitching moment, and z-ordinate at the wing-root trailing edge. 

The total pitching -moment coefficient resulting from n wing-loading distributions 
is given by 

i=n 

= ^ C m iAj (4) 

i=l 

where C mj i denotes the pitching -moment coefficient of the ith loading and is the 
load strength factor for the ith loading. Similarly, the z-ordinate at the wing-root trail- 
ing edge resulting from n wing -loading distributions is given by 
i=n 

z r = £ z r,i A i < 5 > 

i=l 

where z • denotes the z-ordinate at the wing-root trailing edge on the camber surface 
required to support the ith loading. 

If, in addition to a given lift, the constraints of zero pitching moment (C m = 0) and 
given z r are imposed on -the drag -minimization problem, the method of Lagrange multi - 
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pliers yields the following set of equations which establishes the relative strength of each 
loading (see ref. 8): 


x L^L, 1 + X m c m,l 
x L C L, 2 + X m C m,2 


i=n 

+ x z z r,l + Yi C D,li A i = 0 
i=l 

i=n 

+ X z z r 2 + ^ C D,2i A i = 0 
i=l 


i=n 


*L C L,n + X m C m,n + ^z^n + £ C D,ni A i = 0 

i=l 


i=n 


^ C L,i A i ‘ C L,des = 0 


i=l 

i=n 


^ Cm,i A i - 0 


i=l 

i=n 


^ z r,i A i ' 

i=l 


z r = 0 


J 


( 6 ) 


The first n 
k=a 

F = I x k^k 

k=l 


equations of equations (6) are of the form = 0 where 
+ Cj). The expressions <p k = 0 are the constraining relations for lift, 


pitching moment, and z -ordinate at the wing -root trailing edge. These constraining rela- 
tions form the last three equations of equations (6). The load strength factors A^, which 
determine the relative strength of each loading for minimum drag, are solved. for, and the 
camber surface for the minimum-drag combination is determined from 


i=n 

z c (x,y) = ^ z Cji (x,y)Ai 
i=l 


(7) 
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Program Description 

The numerical design method has been programed for the CDC 6600 computer. 

See appendix A for a description of the input required and appendix B for a description 
of the output. This program requires a storage of 55K octal. 

Several options are provided for the user. Any or all of the eight specified load- 
ings shown in figure 5 can be used; also, the pitching -moment constraint (C m = 0) and/or 
the z constraint (specified value of z r ) can be requested along with the lift constraint 
(specified Cj_, )( jes) in the computer program input (see appendix A). 

ILLUSTRATIVE EXAMPLES 
Wing Planforms 

To illustrate and establish the validity of the numerical method, a series of example 
cases is presented for two wing planforms. The results for planform A are presented in 
figures 6 to 11, and the results for planform B are presented in figures 12 to 17. Wing 
planform A is representative of a high Mach cruise vehicle, and for the purpose of demon- 
stration the wing design procedure is applied at M = 3.5; wing planform B is representa- 
tive of a highly maneuverable low -supersonic aircraft, and the wing design procedure is 
applied at M = 1.4. 


Planform A 


Example cases .- Each application presented for wing planform A is denoted by a 
case number, as described in the following table: 


Case 

1 

2 

3 

4 


Description 

Eight loadings, 11 spanwise stations, design for given lift with 

Eight loadings, 27 spanwise stations, design for given lift with 

Eight loadings, 51 spanwise stations, design for given lift with 

Eight loadings, 27 spanwise stations, design for given lift with 

z r constraints 

Three loadings, 27 spanwise stations, design for 
z r constraints 

Case 3 with modified camber at = 0 

b/2 

Eight loadings , 27 spanwise stations , no C m 


given lift with 


Cm 

constraint 

Cm 

constraint 

Cm 

constraint 

Cm 

and 

Cm 

and 


Moment reference location .- Wing planform A is shown in figure 6 with the moment 
reference location which was selected for all seven cases. This moment reference loca- 
tion represents the center of gravity that will provide for a subsonic static margin of 
0.05c as calculated by the Langley subsonic aerodynamic program A2794 (see ref. 9). 
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Program A2794 has given good correlation with experimental static -margin data (see 
ref. 10) and should be a valuable design tool in conjunction with supersonic camber design 
with moment constraint. 

Effect of number of spanwise computation stations .- It was found that the number of 
spanwise computation stations (see JBYMAX in appendix A) significantly affects the solu- 
tion. This is believed to be due to the numerical technique currently used in the program 
to integrate the section lift, drag, and pitching -moment coefficients to obtain the wing lift, 
drag interference, and pitching -moment coefficients ^D,ij> Cniji)* 

Figure 7 illustrates the effect of the number of spanwise computation stations on 
some representative lifting pressure distributions and camber surfaces from the design 
program, and the longitudinal aerodynamic characteristics from the evaluation program 
(see ref. 11). The evaluation program was used to substantiate the results of the design 
program for a Cl d es of 0.10. (Preliminary experimental data indicate that a 
Cl des 0*10 is too high for Mach 3.5 and that theoretically predicted pressures will 
not be realized.) The evaluation program also complements the design program in that 
it provides the drag and pitching-moment characteristics throughout the lift -coefficient 
range for a given camber surface. Figure 7(c) indicates that the moment constraint pro- 
duced zero pitching -moment coefficient at C L des * or each case shown, and the drag 
decreased slightly with increased number of spanwise computation stations. 

Although the parameters have been nondimensionalized, the camber lines at the top 

of figures 7(b), 9(b), and 10(a) [ variation of ■ Z £ with x/l ) are plotted to such a 

V zc L,des / 

scale that the slopes represent the true camber slopes for C^des = 0.10. The camber- 
line slopes for 11 spanwise computation stations are not as severe as those for 27 and 
51 stations, but the results shown in figure 8 indicate that at least 27 stations should be 
used. Figure 8 shows the variation of drag -due -to -lift factor with number of spanwise 
computation stations. Results from the design program show no significant variation in 
drag -due -to -lift factor between 27 and 51 spanwise computation stations. Also, the agree- 
ment between the design and evaluation programs is better for 27 and 51 stations than for 
11 stations. It is believed that the differences between the two programs are accounted 
for by the fact that the evaluation program must extract lifting -element slopes in the 
streamwise direction from the camber-surface y,z ordinates which were provided by 
the design program. These slopes are not necessarily in exact agreement with those 
used in the design program. 

Use of z r constraint and effect of number of loadings used.- Since the camber 
surfaces shown in figure 7(b) for cases 2 and 3 drop excessively at the root chord, a z r 
constraint value of -1.27ZCl des was ac *ded for cases 4 and 5. Figure 9 illustrated the 
effect of the number of specified loadings selected when the moment and z r constraints 
are applied. For case 4 all eight loadings were used, with the specified area for the 
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eighth loading defined by the shaded area in figure 2. For case 5 the three loadings used 
in the method of reference 3 were applied: uniform, linear chordwise, and linear span- 
wise. Figure 9(b) indicates that the solution provided the requested z r value for both 
cases, and the camber lines for case 4 are more reasonable than those for either case 5 
or case 3. Figure 9(c) indicates that the drag for case 4 is substantially lower than that 
for case 5 and only slightly higher than that for case 3 (no zr constraint). 


Effect of modifying the minimum -drag root camber . - Since the root camber slopes 

for case 4 still are somewhat severe, the camber line at = 0.02 was substituted at 

b/2 

= 0. The resulting design, designated case 6, is compared with the results for the 
unconstrained design (case 7) in figure 10. Case 6 is clearly superior to case 7 because 
the camber slopes are much less severe, pitching moment is zero at the design Cl 
( see fig. 10(b)), and the drag is only slightly higher. Therefore it appears that engineer- 
ing judgment is justifiable in the wing design process. The modification method of sub- 
stituting the = 0.02 camber line at = 0 is purely arbitrary and may not be the 
best method in this particular case. 


Figure 11 summarizes the drag -due -to -lift factors for cases 3 to 7. Note that 
C-^/c L des * rom the evaluation program is lower for case 6 than for case 4, indicating 
that the modification to the root chord was not detrimental. 


Planform B 

Wing planform B, shown in figure 12, represents a planform that might be employed 
on a highly maneuverable M = 1.4 airplane. The results of the design method applied 
at M = 1.4 to this planform serve to complement those obtained at M = 3.5 for plan- 
form A. 

Example cases .- For planform B the design procedure was applied with and without 
the pitching -moment constraint for loadings 1 to 3 and loadings 1 to 7. Each application 
presented for wing planform B is denoted by a case number, as described in the following 
table: 

Case Description 

1 Three loadings, design for given lift • 

2 Seven loadings, design for given lift 

3 Three loadings, design for given lift and zero pitching moment 

4 Seven loadings, design for given lift and zero pitching moment 

The eighth loading and the z r -constraint were not employed in this example, in order to 
present an application of the design method which differed from that used for planform A. 
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Moment reference location .- In order to exercise the pitching -moment constraint 
in a meaningful manner, the moment-center location of 0.662Z was determined by per- 
forming a weight -and-balance analysis on a complete airplane configuration employing 
the wing planform of this example. 


Effect of spanwise computation stations .- Figure 13 shows the variation of the drag- 
due-to-lift factor Cj^/c^ ^ gg with number of spanwise computation stations. The cir- 
cular symbols correspond to values obtained from the design procedure of this report; 
the square symbols were obtained by applying the evaluation procedure of reference 11 to 
the minimum -drag camber surface at the design lift coefficient. Results from the design 
program show no significant variation of drag -due -to -lift factor with number of spanwise 
computation stations; however, design -program and evaluation -program results are in 
best agreement for the maximum number of spanwise computation stations. The remain- 
ing discussion is predicated upon results obtained with 51 spanwise computation stations. 


Effect of number of specified loadings and moment constraint . - The effect of num- 
ber of loading distributions on the optimum lifting pressure distribution and the camber 
surface for a wing designed for given lift and no moment constraint is shown in figure 14. 
The lifting pressure distributions which correspond to the optimum combination of load- 
ings are shown in figure 14(a). The paramount differences between the pressure distribu- 
tions are in the complexity and the pressure gradients in the chordwise direction. The 
maximum pressure coefficient shown for seven loadings is over twice that obtained for 
three loadings; however, it is not known whether this type of loading distribution is practi- 
cal. In order to assess the practicality of the more complex loading distribution it is 
necessary to select a particular wing thickness distribution and then consider the total 
combination of thickness and lifting pressures. Such consideration is beyond the scope 
of this report, and only lifting pressure distributions are presented here. 


The camber surfaces shown in figure 14(b) have similar shapes except at tX- = 0. 

u/ 2i 

In this region the more complex loading distribution produces a camber surface which is 
also more complex but possesses a smaller overall z c variation from leading to trail- 


ing edge. For conventional aircraft 


y . .. 
b72" 


0 is in the fuselage, and it is usually difficult 


to achieve substantial amounts of effective camber there; thus, the root -chord camber 
produced by seven loadings is more practical from this consideration. 


For the wing designed for given lift and zero pitching moment, lifting pressure dis- 
tributions and camber surfaces are shown in figure 15. A comparison of the camber sur- 
faces in figure 15(b) again shows that the one corresponding to seven loadings has the 
smaller z c difference between the leading and trailing edges of the root chord, but the 
camber-surface slopes dz c /dx are much greater near the leading edge for seven load- 
ings than for three loadings. These differences are reflected in the pressure distribu- 
tions shown in figure 15(a). 
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For design lift coefficients of 0, 0.1, 0.2, and 0.3, minimum -drag wing camber 
surfaces were generated at M = 1.4 with and without moment constraints, employing 
three loadings and seven loadings. Each of these camber surfaces was evaluated by the 
method of reference 11, and the resulting aerodynamic coefficients of pitching moment 
and of drag due to lift are presented in figure 16 as functions of lift coefficient. 

Aerodynamic coefficients for wing camber surfaces designed for given lift by using 
three loadings (case 1) and seven loadings (case 2) are presented in figures 16(a) and 
16(b), respectively. The drag at the design lift is slightly lower for case 2 than for case 1 
however, the pitching-moment characteristics are substantially different. At design lift, 
three loadings produced a wing camber surface which had a relatively large negative 
pitching moment. Trimming this moment out with a control surface would result in a 
trim drag increment in addition to the drag shown in figure 16(a). For this particular 
planform, the use of seven loadings without moment constraint essentially produced a 
self -trimming wing with a small positive moment at the design lift (unlike case 7 for 
planform A). Figures 16(c) and 16(d) show aerodynamic characteristics of wings designed 
for zero moment at given lift by wing three loadings (case 3) and seven loadings (case 4), 
respectively. All the twisted and cambered wings of figure 16 had better aerodynamic 
characteristics than the flat wing (i.e., design lift of zero) at the higher lift coefficients. 
Although the wings were designed for zero moment, application of the evaluation proce- 
dure resulted in moments near zero but not exactly zero. As previously explained, this 
discrepancy between the wing design and wing evaluation procedures is attributed to the 
numerical techniques employed in evaluating the camber -surface slopes. The lift -drag 
characteristics shown in figure 16(d) for case 4 are essentially the same as those shown 
in figure 16(b) for case 2 and in figure 16(c) for case 3. 

The pitching -moment curves for case 2 are displaced somewhat in the positive 
direction as compared with the curves for the other cases in figure 16. In contrast to 
example A, the application of seven loadings with no pitching-moment constraint (case 2) 
produced an optimum wing possessing a small positive moment at the design lift coeffi- 
cient. In this situation, constraining the pitching moment to be zero would not provide 
significant improvements in the moment characteristics and might even be detrimental to 
the lift -drag characteristics. These two examples clearly illustrate that the benefits to 
be gained by use of v the various options provided with the design program are dependent 
on the planform. •' - - • =■ 

The drag -due -to -lift factor, which serves as a measure of the wing efficiency, is 
summarized in figure 17 for the four cases of planform B. In all cases the drag-due-to- 
lift factors obtained from the evaluation program are larger than those obtained from the 
design program. The factors are smaller for seven loadings than for three loadings, and 
for three loadings the addition of the moment constraint results in an increase in drag- 
due-to-lift factor, but for seven loadings the drag -due -to-lift factors are the same. 
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Possible Extensions 


Although the present design method has been applied only to a wing alone, it should 
be applicable to multiple lifting surfaces such as wing- tail and wing- canard combinations. 
The evaluation capability for wing-tail and wing- canard combinations exists and has been 
successfully demonstrated (ref. 11); similar success is anticipated for the design of these 
combinations. The significance of such a design capability can be appreciated by noting 
that the usual procedure is to design an optimum wing and then add a horizontal tail in ' 
order to trim the configuration. However, if a separate control surface is to be used in 
combination with a wing, the most direct procedure for arriving at an optimum combina- 
tion would be to incorporate the control- surface effects into the initial wing design process. 

This approach solves the complex problem of determining how much trim moment 
should be obtained from wing twist and camber and how much from the control surface. 

This type of solution also determines whether the lift on the control surface should be 
positive or negative for a given center-of-gravity position. 

CONCLUSIONS 

A numerical method based on linearized theory has been developed and programed 
for the CDC 6600 computer to design a minimum-drag supersonic-wing camber surface 
of arbitrary planform for a given lift, with options for constraining the pitching moment 
and/or the surface deformation at the trailing edge of the root chord and for selecting 
any desired combination of eight specified wing-loading distributions to be employed in 
the optimization procedure. The application of this method to two examples led to the 
following conclusions: 

1. For the two wing planforms investigated, the drag penalties for designing wings 
to be self -trimming and to have a reasonable camber surface are relatively small. 

2. The use of seven and eight loadings consistently produced camber surfaces having 
lower drag than camber surfaces produced by three loadings. 

3. The benefits obtained by using constraints on the pitching moment and the root- 
chord trailing-edge ordinate are dependent on the planform. 

4. The number of spanwise computation stations required for a valid solution is 
dependent on the number of specified loadings used. 

5. Replacing the wing-root camber line computed by the program with a camber line 
from the adjacent outboard span station does not significantly alter the drag and can pro- 
vide a more reasonable camber surface. 

Langley Research Center, - ' ; ■ j 

National Aeronautics and Space Administration, 

Hampton, Va., November 20, 1972. 
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APPENDIX A 


INPUT DATA 

The purpose of this appendix is to provide a description of the input required for the 
computer program. All input is in NAMELIST format (ref. 12). A description of the 
required inputs and the FORTRAN variables used by the computer program are as follows: 


FORTRAN 

Variable 


Description 

One card of Hollerith identifying information 


$ INPUT 


Arbitrary name required by the loading routine to define first 
input -data block 


XM 


Free -stream Mach number 


XMAX 

JBYMAX* 


Largest value of x in wing definition 

Total number of spanwise stations for camber-surface 
computation (JBYMAX S 51) 


JBYS(l) 


Table of spanwise stations for camber-surface computation, 
integer values (0 i JBYS = NON) 


NON 3 * 


Number of semispan grid elements selected to represent wing 
(NON g 100) 


XLEO 

XTEO 

RATIO 


x- coordinate of wing leading edge at y = 0 
x-coordinate of wing trailing edge at y = 0 

Factor used to proportion z c -ordinates for given design lift 

coefficient C Mes^>) 

\ NON / 


NOM* 


Number of Mach numbers other than XM 


REFAR 


Reference area of wing 
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FORTRAN 

Variable 

SPAN 

ICG* 

XCG 

IZREST* 

ZMIN 

LOAD1* 

LOAD2* 

LOAD3* 

LOAD4* 

LOAD 5* 

LOAD 6* 

LOAD 7* 

LOAD8* 


APPENDIX A - Continued 


Description 


Total span of wing 

Selection for moment constraint (ICG = 1 to exercise moment 
constraint, input XCG) 

x-distance from XLEO to pitching -moment center 

Selector for zr constraint (IZREST = 1 to exercise z r con- 
straint, input ZMIN) 

z -distance from horizontal reference plane to trailing edge of 
root chord (z r ) 

Selector for loading 1 (LOAD1 = 1 applies uniform loading) 

Selector for loading 2 (LOAD2 = 1 applies linear chordwise 
loading) 

Selector for loading 3 (LOAD3 = 1 applies linear spanwise 
loading) 

Selector for loading 4 (LOAD4 = 1 applies quadratic spanwise 
loading) 

Selector for loading 5 (LOAD5 = 1 applies parabolic chordwise 
loading) 

Selector for loading 6 (LOAD6 = 1 applies quadratic chordwise 
loading) 

Selector for loading 7 (LOAD7 = 1 applies cubic chordwise 
loading) 

Selector for loading 8 (LOAD8 = 1 applies a linear chordwise 
loading on the specified area) 
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FORTRAN 

Variable 

APUNCH 


APLOT 

TMACH(l) 

TXLE(l) 

TXTE(l) 

TYPEX 

NLEX* 

NTEX* 

TBLEX(l) 

TBLEY(l) 

TBTEX(l) 

TBTEY(l) 

$ 


APPENDIX A - Continued 

Description 


Variable to control punch output. APUNCH = 0. for no punch output, 
APUNCH = 1. to punch, z c -ordinates in format (4(E18.8,1H)) com- 
patible with input to program D2340, APUNCH = 2. to punch 
z c -ordinates in format (10(F7.3),2X,5HTZORD,I2) compatible 
with input to D2500 wave-drag program, APUNCH = 3. to punch 
z c -ordinates in both of the above formats. 

Variable to control plotting of wing camber. APLOT = 0. for no 
plots, APLOT =-l. to request plots. ; 

Table of other Mach numbers corresponding to NOM 

Table of x-coordinates of wing leading edge at successive span 
stations of y = where N = 1 to NON; used only when 

TYPEX = O'. 

Table of x-coordinates of wing trailing edge at successive span 

stations of y = where N = 1 to NON; used only when 

NON 

TYPEX = 0. 

If TYPEX = 0. input TXLE and TXTE tables and omit remaining 
input variables in this data block; if TYPEX = 1. TXLE and TXTE 
tables will be generated from the following input information 

Number of leading-edge definition points (=50) 

Number of trailing -edge definition points (=50) 

Table of NLEX wing leading -edge x-coordinates 

Table of : NLEX wing le ading -edge y- coordinates 

Table of NTEX wing trailing-edge x-coordinates 

Table of NTEX. wing trailing-edge y-coordinates 

Denotes end of input-data block 
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APPENDIX A - Continued 


FORTRAN 

Variable 


The following input -data block is required when applying loading 8 (i.e., LOAD 8 = 1). 

Description 


$INPT8 


Arbitrary name required by the loading routine to define second 
input -data block 


XLEOLOD 

TXLOAD(l) 


x-coordinate of specified-area leading edge at y = 0 

Table of x-coordinates of specified-area leading edge at success- 
ive span stations of y = where N = 1 to NON; used only 

when TYPEX = 0. 


NX8* 


TBX8(1) 

$ 


Number of leading-edge definition points describing specified area 
(=50); TYPEX = 1. 

Table of NX8 x-coordinates at specified-area leading edge 
Denotes end of input-data block 
Fixed -point number 
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APPENDIX B 


OUTPUT DATA 

The input data, in real dimensions, are printed out. All other dimensional output 

numbers are in program dimensions (multiplied by NON//3(b/2)) except those multiplied 

by RATIO. The input data in program dimensions are next printed out along with the 

pressure -coefficient terms for the eight loadings: 

CPU Lifting pressure coefficient of uniform load =1.0 

DCPDBY Lifting pressure coefficient of linear spanwise loading per unit j3y = 0.04 

DCPDXP Lifting pressure coefficient of linear chordwise loading 

per unit x r = 0.04 

DCPDBY2 lifting pressure coefficient of parabolic spanwise loading 
per unit /3 2 y 2 = 0.04 

DCPDX5 Lifting pressure coefficient of parabolic chordwise loading 

per unit (x') 2 = 0.0004 

DCPDX6 Lifting pressure coefficient of quadratic chordwise loading 

per unit (x')(x T - c) = 0.0016 

DCPDX7 Lifting pressure coefficient of cubic chordwise loading 

per unit (x') 2 (1.5c - x') = 0.00008 

DCPDX8 Lifting pressure coefficient of linear chordwise loading per unit x' 

on the specified area = 0.04 


The drag coefficients of interference between each loading on the camber surface 
and the other loadings are next printed, along with the following values for each loading: 


AREA 

CBAR 


pb/2 

Right-hand wing area, \ c dy 

J 0 


Wing mean aerodynamic chord, 


f b/2 c 2 dy/AREA 
J 0 
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APPENDIX B - Continued 
r b/2 

CD Wing drag coefficient, \ (CDS)c dy/AREA 

J 0 

pb/2 

CL Wing lift coefficient, j (CLS)c dy/AREA 

CMCG Wing pitching -moment coefficient about XCG, 

f b / 2 (CMS)c 2 dy/^REA(CBAR) 

The force coefficients of the specified and interference loadings are printed with 
other terms related to the solution of the simultaneous equations: 

CDij Drag coefficient of interference between loadings i and j 

(see eq. (20) of ref. 3) 

CLi Lift coefficient of ith loading 

DETERM Value of the determinant of the coefficient matrix used in solving 

equations (3) 

Ai Proportionality factor for ith loading to produce minimum -drag combina 

tion having total lift coefficient of 1.0 

LAMBDA1 of simultaneous equations ^C ^ constraint) 

LAMBDA2 \ m of simultaneous equations (moment constraint) 

LAMBDA3 A z of simultaneous equations (z constraint) 

The wing data for the minimum -drag loading are printed for each semispan station 
selected (see JBYS(l) in appendix A). The calculated ZRATIO values are interpolated at 
standard NACA chordwise stations. 

YB2 Semispan station 

XRATIO x-distance to the center of the grid element times RATIO 

L x-distance to trailing edge of grid element 
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APPENDIX B - Concluded 

DZDX Streamwise camber -surface slope at x or x' , 9z c /9x 

ZRATIO z c times RATIO, value of camber surface ordinate at x or x' , 

RATI <(^' 

LSTAR Grid element number, ascending order in x-direction 

CP/CL Lifting pressure coefficient at trailing edge of grid element 

divided by wing lift coefficient 

XPOC x* expressed as fraction of local chord, x’/CHORD 


CHORD 

Chord length at selected YB2, (TXTE - TXLE) JByg 

r 

= c 
1 

9z„ /„*\ 

CDS 

Section axial -force coefficient at selected YB2, -\ 

J C 

( AC P 

— — d (— ) 

9x \cj 

CLS 

Section lift coefficient at selected YB2, J AC p d^- 

r) . 


CMS 

Section pitching -moment coefficient about XCG, -f 

1 

AC p 

x - XCG , 
c a 


J 0 


The drag -due -to -lift factor (CD/CL**2) based on REFAR, and the values of AREA, 
CBAR, CD, CL, and CMCG for the minimum -drag wing, are printed last. 
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Figure 1.- Potential use of wing twist and camber to reduce trim drag. 
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Specified region 









Figure 8.- Variation of drag -due -to -lift factor with number of spanwise computation 

stations for planform A wing. 
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(a) Lifting pressure distributions. 

Figure 9.- Effect of number of specified loadings on design of planform A wing with 
moment and z r constraint. z r = -1.27ZCLd es ; M = 3.50. 
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Case 5, Case 4, Case 6, Case 3, . Case 7, 

moment moment modified without z r without 

and z r and z r case 4 constraint moment 

constraint constraint and z r 

constraint 

Figure 11.- Summary of effect of design program options on drag -due -to -lift factor 

for planform A wings. 


36 




c Design program 
□ Evaluatio n p rog ra m 


Case 1 
3 loadings 



Case 2 
7 loadings 

. 25 r n No moment constraint 



Case 3 
3 loadings 


• 25r Moment constraint 



10 20 30 40 50 

Spanwise computation stations 


Figure 13.- Variation of drag -due -to -lift factor with number of spanwise 
computation stations for planform B wings. 
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(a) Lifting pressure distributions. 

Figure 14.- Effect of number of specified loadings on design of planform B wings 
for given lift without moment constraint. M = 1.4. 






Case 3 


y 3 loadings 




(a) Lifting pressure distributions. 

Figure 15.- Effect of number of specified loadings on design of planform B wings 
for given lift with moment constraint. M = 1.4. 
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